An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2.

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

The importance of accessory protein variants in the pathogenicity of SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS- CoV-2) with an estimated fatality rate of less than 1%. The SARS-CoV-2 accessory proteins ORF3a, ORF6, ORF7a, ORF7b, ORF8, and ORF10 possess putative functions to manipulate host immune mechanisms. These involve interferons, which appear as a consensus function, immune signaling receptor NLRP3 (NLR family pyrin domain-containing 3) inflammasome, and inflammatory cytokines such as interleukin 1ß (IL-1ß) and are critical in COVID-19 pathology. Outspread variations of each of the six accessory proteins were observed across six continents of all complete SARS-CoV-2 proteomes based on the data reported before November 2020. A decreasing order of percentage of unique variations in the accessory proteins was determined as ORF3a > ORF8 > ORF7a > ORF6 > ORF10 > ORF7b across all continents. The highest and lowest unique variations of ORF3a were observed in South America and Oceania, respectively. These findings suggest that the wide variations in accessory proteins seem to affect the pathogenicity of SARS-CoV-2.

Structural and Functional Genomics of the Resistance of Cacao to Phytophthora palmivora.

Black pod disease, caused by Phytophthora spp., is one of the main diseases that attack cocoa plantations. This study validated, by association mapping, 29 SSR molecular markers flanking to QTL (Quantitative Trait Loci) associated with Phytophthora palmivora Butler (Butler) (PP) resistance, in three local ancient varieties of the Bahia (Comum, Pará, and Maranhão), varieties that have a high potential in the production of gourmet chocolate. Four SSR loci associated with resistance to PP were detected, two on chromosome 8, explaining 7.43% and 3.72% of the Phenotypic Variation (%PV), one on chromosome 2 explaining 2.71%PV and one on chromosome 3 explaining 1.93%PV. A functional domains-based annotation was carried out, in two Theobroma cacao (CRIOLLO and MATINA) reference genomes, of 20 QTL regions associated with cocoa resistance to the pathogen. It was identified 164 (genome CRIOLLO) and 160 (genome MATINA) candidate genes, hypothetically involved in the recognition and activation of responses in the interaction with the pathogen. Genomic regions rich in genes with Coiled-coils (CC), nucleotide binding sites (NBS) and Leucine-rich repeat (LRR) domains were identified on chromosomes 1, 3, 6, 8, and 10, likewise, regions rich in Receptor-like Kinase domain (RLK) and Ginkbilobin2 (GNK2) domains were identified in chromosomes 4 and 6.

RRGPredictor, a set-theory-based tool for predicting pathogen-associated molecular pattern receptors (PRRs) and resistance (R) proteins from plants.

In plant-pathogen interactions, plant immunity through pathogen-associated molecular pattern receptors (PAMPs) and R proteins, also called pattern recognition receptors (PRRs), occurs in different ways depending on both plant and pathogen species. The use and search for a structural pattern based on the presence and absence of characteristic domains, regardless of their disposition within a sequence, could be efficient in identifying PRRs proteins. Here, we develop a method mainly based on text mining and set theory to identify PRR and R genes that classify them into 13 categories based on the presence and absence of the main domains. Analyzing 24 plant and algae genomes, we showed that the RRGPredictor was more efficient, specific and sensitive than other tools already available, and identified PRR proteins with variations in size and in domain distribution throughout the sequence. Besides an easy identification of new plant PRRs proteins, RRGPredictor provided a low computational cost.

Involvement of structurally distinct cupuassu chitinases and osmotin in plant resistance to the fungus Moniliophthora perniciosa.

The cupuassu tree (Theobroma grandiflorum) is a crop of great economic importance to Brazil, mainly for its pulp and seeds, which are used in food industry. However, cupuassu fruit production is threatened by witches' broom disease caused by the fungus Moniliophthora perniciosa. As elements of its defense mechanisms, the plant can produce and accumulate pathogenesis-related (PR) proteins such as chitinases and osmotins. Here, we identified three cupuassu PR proteins (TgPR3, TgPR5 and TgPR8) from cupuassu-M. perniciosa interaction RNA-seq data. TgPR3 and TgPR8 corresponded to chitinases, and TgPR5 to osmotin; they are phylogenetically related to cacao and to Arabidopsis PR sequences involved in biotic and abiotic stress. The TgPR proteins' tridimensional structure was obtained through homology modeling, and molecular docking with chitin and chitosan showed that the TgPR proteins can interact with both cell wall molecules and presented a higher affinity for chitosan. TgPR gene expression was analyzed by RT-qPCR on resistant and susceptible cupuassu genotypes infected by M. perniciosa at 8, 24, 48 and 72 h after infection (hai). The TgPR genes showed higher expression in resistant plants compared to the susceptible ones, mainly for TgPR5 at 8 and 24 hai, while the expression was lower in the susceptible cupuassu plants. To our knowledge, this is the first in silico and in vitro reports of cupuassu PR protein. The data suggested that TgPRs could be involved in recognizing mechanisms of the plant's innate immune system through chitin receptors. Our results also suggest a putative role of chitinase/chitosanase for the TgPR5/osmotin.

HVA22 from citrus: A small gene family whose some members are involved in plant response to abiotic stress.

The HVA22 gene has been isolated for the first time from the aleurone layer of barley (Hordeum vulgare). Here, we characterized the HVA22 family from citrus (C. clementina and C. sinensis). Twelve genes, 6 in each species, were identified as well as duplication events for some of them. The ORF size ranged from 235 to 804 bp and the protein molecular weight from 94 to 267 kDa. All the citrus HVA22 protein presented transmembrane location and conserved TB2/DP1/HVA22 region. Phylogenetic and gene expression analyses suggested that some citrus HVA22 play a role in flower and fruit development, and that gene expression may be regulated by hormone or environmental conditions. Other regulation levels were also predicted, such as alternative splicing and post-translational modifications. The overall data indicated that citrus HVA22 may be involved in vesicular traffic in stressed cells, and that CcHVA22d could be involved in dehydration tolerance.

Genome-wide identification and characterization of cacao WRKY transcription factors and analysis of their expression in response to witches' broom disease.

Transcriptional regulation, led by transcription factors (TFs) such as those of the WRKY family, is a mechanism used by the organism to enhance or repress gene expression in response to stimuli. Here, we report on the genome-wide analysis of the Theobroma cacao WRKY TF family and also investigate the expression of WRKY genes in cacao infected by the fungus Moniliophthora perniciosa. In the cacao genome, 61 non-redundant WRKY sequences were found and classified in three groups (I to III) according to the WRKY and zinc-finger motif types. The 61 putative WRKY sequences were distributed on the 10 cacao chromosomes and 24 of them came from duplication events. The sequences were phylogenetically organized according to the general WRKY groups. The phylogenetic analysis revealed that subgroups IIa and IIb are sister groups and share a common ancestor, as well as subgroups IId and IIe. The most divergent groups according to the plant origin were IIc and III. According to the phylogenetic analysis, 7 TcWRKY genes were selected and analyzed by RT-qPCR in susceptible and resistant cacao plants infected (or not) with M. perniciosa. Some TcWRKY genes presented interesting responses to M. perniciosa such as Tc01_p014750/Tc06_p013130/AtWRKY28, Tc09_p001530/Tc06_p004420/AtWRKY40, Tc04_p016130/AtWRKY54 and Tc10_p016570/ AtWRKY70. Our results can help to select appropriate candidate genes for further characterization in cacao or in other Theobroma species.

Alternative oxidase (AOX) constitutes a small family of proteins in Citrus clementina and Citrus sinensis L. Osb.

The alternative oxidase (AOX) protein is present in plants, fungi, protozoa and some invertebrates. It is involved in the mitochondrial respiratory chain, providing an alternative route for the transport of electrons, leading to the reduction of oxygen to form water. The present study aimed to characterize the family of AOX genes in mandarin (Citrus clementina) and sweet orange (Citrus sinensis) at nucleotide and protein levels, including promoter analysis, phylogenetic analysis and C. sinensis gene expression. This study also aimed to do the homology modeling of one AOX isoform (CcAOXd). Moreover, the molecular docking of the CcAOXd protein with the ubiquinone (UQ) was performed. Four AOX genes were identified in each citrus species. These genes have an open reading frame (ORF) ranging from 852 bp to 1150 bp and a number of exons ranging from 4 to 9. The 1500 bp-upstream region of each AOX gene contained regulatory cis-elements related to internal and external response factors. CsAOX genes showed a differential expression in citrus tissues. All AOX proteins were predicted to be located in mitochondria. They contained the conserved motifs LET, NERMHL, LEEEA and RADE-H as well as several putative post-translational modification sites. The CcAOXd protein was modeled by homology to the AOX of Trypanosona brucei (45% of identity). The 3-D structure of CcAOXd showed the presence of two hydrophobic helices that could be involved in the anchoring of the protein in the inner mitochondrial membrane. The active site of the protein is located in a hydrophobic environment deep inside the AOX structure and contains a diiron center. The molecular docking of CcAOXd with UQ showed that the binding site is a recessed pocket formed by the helices and submerged in the membrane. These data are important for future functional studies of citrus AOX genes and/or proteins, as well as for biotechnological approaches leading to AOX inhibition using UQ homologs.

Selection of Reference Genes for Expression Study in Pulp and Seeds of Theobroma grandiflorum (Willd. ex Spreng.) Schum.

Cupuassu (Theobroma grandiflorum [Willd. ex Spreng.] Schum) is a species of high economic importance in Brazil with great potential at international level due to the multiple uses of both its seeds and pulp in the industry of sweets and cosmetics. For this reason, the cupuassu breeding program focused on the selection of genotypes with high pulp and seed quality-selection associated with the understanding of the mechanisms involved in fruit formation. Gene expression is one of the most used approaches related to such understanding. In this sense, quantitative real-time PCR (qPCR) is a powerful tool, since it rapidly and reliably quantifies gene expression levels across different experimental conditions. The analysis by qPCR and the correct interpretation of data depend on signal normalization using reference genes, i.e. genes presenting a uniform pattern of expression in the analyzed samples. Here, we selected and analyzed the expression of five genes from cupuassu (ACP, ACT, GAPDH, MDH, TUB) to be used as candidates for reference genes on pulp and seed of young, maturing and mature cupuassu fruits. The evaluation of the gene expression stability was obtained using the NormFinder, geNorm and BestKeeper programs. In general, our results indicated that the GAPDH and MDH genes constituted the best combination as reference genes to analyze the expression of cupuassu samples. To our knowledge, this is the first report of reference gene definition in cupuassu, and these results will support subsequent analysis related to gene expression studies in cupuassu plants subjected to different biotic or abiotic conditions as well as serve as a tool for diversity analysis based on pulp and seed quality.